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Abstract. - We investigate the effect of Dzyaloshinskii-Moriya (DM) interactions on torque 
measurements of quantum magnets with magnetization plateaux in the context of a frustrated 
spin-1/2 ladder. Using extensive DMRG simulations, we show that the DM contribution to the 
torque is peaked at the critical fields, and that the total torque is non-monotonous if the DM 
interaction is large enough compared to the g-tensor anisotropy. More remarkably, if the DM 
vectors point in a principal direction of the g-tensor, torque measurements close to this direction 
will show well defined peaks even for small DM interaction, leading to a very sensitive way to 
detect the critical fields. We propose to test this effect in the two-dimensional plateau system 

SrCu2(B03)2. 
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The investigation of quantum magnets in high mag- 
netic fields is a very active field of research thanks to a 
ilumber of recent and remarkable discoveries ranging from 
Bose-Einstein condensation [1] to magnetization plateaux 
[2,3], and to the on-going search for the analog of su- 
persolid phases [4-7]. A central piece of information is 
provided by the magnetization as a function of the exter- 
nal field. However, technical requirements have lead many 
groups to measure the torque rather than the magnetiza- 
tion. In SU(2) invariant magnets, it does not matter since 
the torque is rigorously proportional to the magnetization 
along the field. In the presence of anisotropic interactions, 
such as Dzyaloshinskii-Moriya (DM) interactions [8], this 
is no longer the case and an additional response is ob- 
tained which needs to be considered carefully. This DM 
contribution to the torque has already been investigated 
in the context of molecular magnets [9,10] and of the spin 
ladder compound Cu(Hp)Cl [11]. In the latter case, the 
DM interactions were assumed not to compete with the 
exchange, so that their additional effect on the torque is a 
relatively smooth contribution in the intermediate phase. 
Recent torque measurements on SrCu2 (603)2, a quasi- 
two-dimensional quantum magnet with several magneti- 
zation plateaux [2,3], call for further investigation of the 
issue. In particular, the field dependence of the torque 
reported in refs. [12, 13] is quite different, especially re- 
garding the behavior inside the plateaux. In Sebastian 




Fig. 1: (Colour on-line) Ladder system under consideration. 
Red arrows indicate DM vectors, the magnetic field lies in the 
{y, 2) plane, but is tilted by an angle a with respect to the 
2- axis. 



et aZ.'s data, the torque is never flat but increases quite 
significantly inside the plateaux, while in Levy et aL's 
data, it actually decreases inside the 1/8 plateau. Since 
DM interactions have been unambiguously identified in 
SrCu2(B03)2 [14-16], the torque is not expected to be 
simply proportional to the magnetization, which leaves the 
door open to anomalous behavior depending on the details 
of the experiments, in particular the orientation of the 
field. However, the field dependence of the torque around 
magnetization plateaux has never been investigated so far. 

In this Letter, we investigate this issue in the context of 
the frustrated ladder depicted in fig. 1 and defined by the 
Hamiltonian {h=\) 
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Fig. 2: (Colour on-line) DMRG results for the uniform magne- 
tization (a) parallel and (b) perpendicular to H for D — 0.03 
and 5g = for a = 30° or 80°, respectively, and — Jx =0.1 
for systems with 25 rungs. In (b), the magnetization is scaled 
by sin a cos a, the angular dependence of the torque of a single 
dimer with DM anisotropy away from the critical field. 



+J\\ ■ Sj+i,i + Sj^2 ■ 5^+1,2 

j 



s 



(1) 
(2) 



This frustrated ladder has a magnetization plateau at 1 /2 
if J|| and Jx are not too different (their ratio should be 
between 1/3 and 3 in the limit Jy, Jx <C J) [17], and it 
is amenable to the density matrix renormalization group 
method (DMRG). It is thus an ideal minimal model to 
study subtle effects related to plateaux in frustrated quan- 
tum magnets. We set J = 1 throughout the paper and we 
will concentrate on the case J|| = Jx for which the plateau 
is largest. We assume the g-tensor g to be diagonal but to 
possess an asymmetry Sg between the x-y and the z com- 
ponent. Due to the DM interaction, the SU(2) symmetry 
of the original Heisenberg model is not present and the 
DMRG calculations can be very demanding, restricting us 
to treat only rather small systems with up to 49 rungs. 

To supplement the finite size DMRG data, we have also 
performed a mean-field calculation which generalizes the 
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Fig. 3: (Colour on-line) Angular dependence of t/H at the 
critical field H = 1 and inside the plateau (H — 1.1) as ob- 
tained from DMRG calculations with L — 25 and L — 49 rungs 
and from the mean-field calculation using the ansatz eq. (3). 



ansatz used in ref. [18] in order to deal with the tilted field. 
It is based on a product of rung wave functions 

[</>), = a\S),+b\T^,),+c\To\),+d\T,), (3) 

with the singlet and triplet eigenstatcs of the dimer on 
rung J given by \S), = l/x/2(| it),- 1 Ti),), l^-i ), = | U 
)j, [To), = l/^/2(| it), + I Ti),), and \T,), = \ tt),- For 
the case J|| ~ Jx and _D = 0, this ansatz provides an exact 
solution [17], and it is expected to provide a good approx- 
imation for the physically relevant case D <^ J . Com- 
parison with the finite system DMRG data shows good 
agreement, giving us confidence that the results presented 
are not affected by strong finite size effects. 

In the following we focus on the magnetic response of 
the system to an external magnetic field with arbitrary 
orientation. First we discuss our results in the presence of 
D only, and afterwards go to the more general case with 
an additional finite 5g. The symmetry analysis for a single 
dimer with DM interaction and 5g = ascertains that the 
results are independent of the absolute orientation of H 
and D as long as the angle between them is the same [19] 
for which reason we confine ourselves to treat only the 
case depicted in fig. 1. If not mentioned otherwise, we 
present results for D = 0.03 and J|| = Jx = 0.1, i.e., 
we restrict ourselves to the case of a strongly frustrated 
ladder, which is relevant when having in mind the frus- 
trated plateau system SrCu2 (603)2. We concentrate on 
the uniform magnetization, which is defined as 
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In general, for finite 5g a magnetization perpendicular to 
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H is induced, causing a torque 

f = ifiu X H, 



(5) 



which gives access to the uniform magnetization perpen- 
dicular to H obtained as = t/H . If D = 0, in order to 
minimize the Zeenian-term in the Hamiltonian, the spins 
align in the direction of the vector iJg. For a better com- 
parison of our results obtained at various angles and for 
different values of i5q, let us introduce the effective field 



(6) 



In the following, when Sg ^ 0, we discuss the dependence 
•mu{H°^). If D = and 6g ^ 0, the components of the 
uniform magnetization parallel and perpendicular to H 
arc obtained as 
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and hence 



Both quantities are proportional to 
and have the same dependence on the magnitude of 
the field so that rriu can be obtained by measuring 

the torque. For a better comparison of the results for dif- 
ferent parameters, we rescale the magnetizations by their 
angular- and g-dependence, i.e., if Sg ^ we analyze 



I ,rescaled 



\/ g"^ + cos^ a 5g{2g + 5g) 



(9) 



^±.rcscalcd 



as a function of 



^Jg"^ +cofi'^ a 5g{2g + 5g) 
' s\nacosa6g{2g + 5g) 



Fig. 4: (Colour on-line) Rescaled uniform magnetization as 
obtained by eqs. (9) and (10) parallel (a) and perpendicular 
(b) to Hior D = 0.03, a = 1°, and 0.005 < 5g < 0.2. The solid 
lines are the results of the mean-field calculation, the points 
the results of DMRG calculations on 49 rungs. The black line 
shows the exact magnetization parallel to H when D = 5g = Q. 
Note that in (a) the mean-field results lie on top of each other. 



which is smoothed out due to the DM interaction, while 
in pronounced peaks appear at the critical fields con- 
fining the plateau. This can be understood by considering 
the angular dependence of the torque T{a). As discussed 
in ref. [19] for the case of a single dimer, it changes dras- 
tically when leaving the limit of perturbatively small H 



= iiBHyJg"^ + cos2 a Sg{2g + Sg). (11) and going to the critical field i?c; one finds 



Before we go to the general case with finite Sg and finite 
D, we first discuss the illustrative case of a system with 
DM interactions only. To do this, it is useful to recall the 
results obtained for an isolated dimer with D ^ 0, Sg = 
[19]. For perturbatively weak fields one finds 
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leading, in general, to a magnetization perpendicular to H. 
Note, however, that in this case it is not possible to deter- 
mine mi by measuring r, since m|(i7) ^ m^{H) due to 
the lack of SU(2) symmetry. This is demonstrated in fig. 
2 where we present DMRG results for the case D = 0.03 
and Sg = 0. As can be seen, the field dependences of m| 
and TO^ differ fundamentally. In mil a plateau is obtained 



This leads to a clear peak in m,^ at He. Note that in fig. 
2(b) we present rescaled by sin a cos a, taking into 
account the angular dependence for <C i?c- As can be 
seen, the same dependence on a is found for all values of H 
sufficiently away from the critical points. In fig. 3 we com- 
pare the angular dependence of at the critical point 
I-IbSHc.i = 1 and inside the plateau for hb^H = 1.1. In 
the latter case, the finite size DMRG and the mean-field 
results are in excellent agreement with each other and with 
the dependence T{a) ~ sin a cos a, and we see that finite 
size effects play a minor role. At the critical field, however, 
we see that also for the ladder systems a completely differ- 
ent angular dependence is obtained. For small a, it follows 
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sin a, while for larger angles m,^ deviates to smaller values. 
Note that at the critical point the results of the DMRG 
and of the mean-field calculations show qualitatively sim- 
ilar behavior, but the mean-field value is smaller than the 
DMRG result. Comparison of DMRG results for 25 and 
49 rungs shows that finite size effects are minimal also in 
this case, so that wc believe that these results should be 
representative for the thermodynamic limit. Despite the 
deviation from the sin a, for large a the resulting at 
the critical field and inside the plateau differ by an order of 
magnitude, providing an explanation for the peaks visible 
in fig. 2(b). 

In the following we test how a finite 5g changes the pic- 
ture obtained for finite D only. In fig. 4 we show our 
results for ladder systems when a = 1° and 0.005 < 5g < 
0.2, i.e. we go to a realistic value of Sg = 10% while keep- 
ing D = 0.03. If 6g < D, the peaks in at the critical 
fields are well visible. For 5g > D, however, the peaks van- 
ish. In this case, the magnetization curves resemble the 
SU(2) symmetric case up to a smoothcning of the plateau, 
and we obtain mi'°"'='''°'^(i/°*'') « m^'^°'''=^i°d(ff'=ff). 

The remarkable angular dependence of the torque ob- 
tained for systems with Sg = has also very important 
consequences in the presence of g-tensor anisotropy. If we 
now keep realistic values of both anisotropics 6g = 10% 
and D = 0.03 and increase a, for angles a > 80° pro- 
nounced peaks in at the critical fields become vis- 
ible again, as can be seen in fig. 5. We find that 

^±,rcscalod(^cff-) _^ ^|l,roscalcd|.^eff ^^IC CaSC with 

DM anisotropy only. This shows that care has to be taken 
when interpreting results of torque measurements, and 
that, in general, it is not possible to safely conclude on 
the field dependence of mi by considering the field depen- 
dence of the torque. 

In conclusion, by investigating the interplay of DM 
and g-tensor anisotropics on the magnetic response of a 
strongly frustrated ladder, we have clarified in which cases 
the torque can be used as a good approximation to the 
magnetization when DM is present. For this to be the case, 
two conditions must be fulfilled: i) The ratio of the DM in- 
teraction to the exchange coupfing D/J should be smaller 
than (or at most comparable to) the g-tensor anisotropy; 
ii) The angle between the magnetic field and the DM vec- 
tor should not be too small. While the first condition is 
often fulfilled (the DM interaction is typically a few per- 
cent of the exchange while the g-tensor anisotropy is about 
10% for Cu^"*"), the second condition may or may not be 
satisfied depending on the experimental conditions. Inter- 
estingly enough, the fact that the field dependence of the 
torque is quite different from that of the magnetization 
when the field is almost parallel to the DM vector could 
be an advantage. Indeed, in this geometry we predict the 
torque to have well pronounced peaks at the critical fields 
that delimitate the plateau. This effect could thus be used 
to locate with a high precision the critical fields, a diffi- 
cult task if only magnetization data are available since 




Fig. 5: (Colour on-line) Same as in fig. 4 when varying the 
angle 1° < Q < 89° for D — 0.03 in the presence of large 
g-tensor anisotropics 5g = 0.2. 

DM interactions can lead to a significant rounding at the 
boundaries of the plateaux. In view of the controversies 
regarding the plateau structure of SrCu2 (603)2, it would 
thus be particularly interesting to perform torque mea- 
surements with the field perpendicular to the layers and 
parallel (or almost parallel) to the intra-dimcr DM vector 
of one type of dimcrs, and it is our hope that the present 
paper will encourage such an investigation. 

* * * 

We acknowledge useful discussions with C. Bcrthier, 
J. Dorier, M. Horvatic, S. Miyahara, K. Penc, I. 
Rousochatzakis, K. P. Schmidt, I. Sheikin, M. Takigawa, 
and T. A. Toth. This work was supported by the Swiss 
National Fund and by MaNEP. 

REFERENCES 

[1] Ch. Riiegg, N. Cavadini, A. Furrer, H.-U. Giidel, 
K. Kramer, H. Mutka, A. Wildes, K. Habicht, and 
R Vorderwisch, Nature 423, 62 - 65 (2003). 

[2] H. Kageyama, K. Yoshimura, R. Stern, N.V. Mushnikov, 
K. Onizuka, M. Kato, K. Kosuge, CP. Slicliter, T. Goto, 
and Y. Ueda, Phys. Rev. Lett. 82, 3168 (1999). 

[3] K. Kodama, M. Takigawa, M. Horvatic, C. Berthier, 
H. Kageyama, Y. Ueda, S. Miyahara, F. Becca, and 
F. Mila, Science 298, 395 (2002). 



Torque anomalies at magnetization plateaux in quantum magnets with Dzyaloshinskii-Moriya interactions 



[4] K.-K. Ng and T.K. Lee, Phys. Rev. Lett. 97, 127204 

(2006) . 

[5] P. Sengupta and CD. Batista, Phys. Rev. Lett. 98, 227201 

(2007) . 

[6] N. Laflorencie and F. Mila, Phys. Rev. Lett. 99, 027202 

(2007) . 

[7] K.P. Schmidt, J. Dorier, A.M. LauchU, and F. Mila, Phys. 

Rev. Lett. 100, 090401 (2008). 
[8] I.E. Dzyaloshinskh, J. Phys. Chem. Solids 4, 241 (1958); 

T. Moriya, Phys. Rev. Lett. 4, 228 (1960). 
[9] F. Cinti, M. Aflronte, and A.G.M. Jansen, Eur. Phys. J. 

B 30, 461-468 (2002). 
[10] V. Lante, L Rousochatzakis, K. Penc, O. Waldmann, and 

F. Mila, arXiv:0810.3837 (2008). 
[11] M. Clemancey, H. Mayaffre, C. Berthier, M. Horvatic, 

J.-B. Fouet, S. Miyahara, F. Mila, B. Chiari, and O. Pi- 

ovesana, Phys. Rev. Lett. 97, 167204 (2006). 
[12] S.E. Sebastian, N. Harrison, P. Sengupta, CD. Batista, 

S. Francoual, E. Palm, T. Murphy, H.A. Dabkowska, and 

B. D. Gaulin, Proc. Natl. Acad. Sci. U.S.A 105, 20157 

(2008) . 

[13] F. Levy, I. Sheikin, C. Berthier, M. Horvatic, M. Taki- 
gawa, H. Kageyama, T. Waki, and Y. Ueda, Europhys. 
Lett. 81, 67004 (2008). 

[14] O. Cepas, K. Kakurai, L.P. Regnault, T. Ziman, 
J. P. Boucher, N. Aso, M. Nishi, H. Kageyama, and 
Y. Ueda, Phys. Rev. Lett. 87, 167205 (2001). 

[15] A. Zorko, D. Arcon, H. van Tol, L.C Brunei, and 
H. Kageyama, Phys. Rev. B 69, 174420 (2004). 

[16] K. Kodama, S. Miyahara, M. Takigawa, M. Horvatic, 

C. Berthier, F. Mila, H. Kageyama, and Y. Ueda, J. Phys.: 
Condens. Matter 17, L61-L68 (2005). 

[17] F. Mila, Eur. Phys. J. B 6, 201 (1998). 

[18] K. Penc, J.-B. Fouet, S. Miyahara, O. Tchernyshyov, and 

F. Mila, Phys. Rev. Lett. 99, 117201 (2007). 
[19] S. Miyahara, J.-B. Fouet, S.R. Manmana, R.M. Noack, 

H. MayafTre, L Sheikin, C. Berthier, and F. Mila, Phys. 

Rev. B 75, 184402 (2007). 



